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ABSTRACT 

X WC-130  aircraft  was  used  as  an  airborne  platform  for  broadband  j 

electric  and  magnetic  field  sensors  to  measure  the  surface  fields  on 
the  airframe  from  lightning  strikes  in  ti.c  7 - 35  km  ranee.  The  wave- 
forms were  recorded  digitally,  with  a sample  interval  of  20  nsec  and  a 

time  window  of  164  microseconds.  Since  the  aircraft  was  above  15,000  ’ 

feet,  and  the  lightning  source  was  relatively  nearby,  the  effects  of 
ground  propagation  were  minimized.  Spectra  have  been  obtained  for 

the  stepped  leader,  first  return  stroke,  and  subsequent  return  stroke  1 

over  the  frequency  range  100  kHz  to  20  MHz.  These  spectra  are 
compared  to  published  data  of  lightning  field  spectra  obtained  from 
ground-based  measurements.  With  airframe  resonances  and  field  enhance- 
ment effects  removed,  the  data  have  two  primary  applications:  , 

(l)  To  add  to  the  data  base  for  the  nearby  lightning  strike 
threat,'  ir\d 

.2 < To  infer  characteristics  of  the  stepped  leader,  first 
return  stroke,  and  subsequent  return  stroke  processes. 

Recommendations  for  these  applications  are  presented. 
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ATTENTION  HAS  BEEN  GIVEN  in  recent  years  to  char- 
acter! zation  of  unwanted  sources  of  electromag- 
netic interference,  such  as  lightning  or  nuclear 
EMP,  that  may  interact  with  aircraft  and  impair 
mission  capability.  The  use  of  non-metallic  air- 
frame materials  and  solid-state  avionics  aggra- 
vates the  problem  by  reducing  electromagnetic 
shielding  and  lowering  damage  or  upset  thresholds 
on  new  generation  aircraft.  Measurements  of 
lightning  fields  and  currents,  indicating  the 
level  of  the  hazards,  have  been  made  on  the 
ground  for  several  years.  Recently,  wide  band- 
width measurements  of  lightning  fields  and  cur- 
rents have  been  recorded  on  aircraft  in  flight  to 
obtain  a more  pertinent  threat  characterization. 
As  a part  of  this  characterization,  it  is  neces- 
sary to  consider  the  electromagnetic  effect  of 
the  aircraft  on  the  measurements  so  that  a gener- 
ic threat  definition  may  be  derived. 

In  1979  the  Air  Force  Flight  Dynamics  Labor- 
atory, in  conjunction  with  the  National  Oceanic 
and  Atmospheric  Administration  (NOAA) , initiated 
a three-year  program  to  measure  the  properties  of 
lightning  encountered  in  flight  (1).  A NOAA  WC- 
130  "hurricane  Hunter"  aircraft  was  fitted  with 
sensors  for  measuring  nearby  and  direct  lightning 
strikes.  A sophisticated  data  acquisition  system 
designed  for  repetitive  wide  bandwidth  transient 
measurements  was  developed  by  the  Micro  Pro  Com- 
pany and  installed  in  the  WC-130.  More  than  one 
hundred  flight  hours  were  logged  over  the  three- 
year  period  and  several  hundred  wide  bandwidth 
lightning  measurements  were  recorded.  Progress 
and  results  of  this  research  are  documented  in 
references  (1)  through  (4). 

In  this  paper,  we  have  analyzed  the  frequen- 
cy characteristics  of  three  lightning  transients 
— stepped  leader  pulses,  first  return  strokes, 
and  subsequent  return  strokes,  recorded  on  the 
WC-130.  Wide  bandwidth  magnetic  field  measure- 
ments were  Fourier  transformed,  averaged,  and 
normalized  for  range.  The  results  are  compared 
with  published  frequency  spectra  from  similar 
field  measurements  obtained  on  the  ground. 

INSTRUMENTATION 

The  sensors  and  instrumentation  used  on  the 
WC-130  have  been  described  in  detail  in  reference 
(4).  The  specific  equipment  used  to  obtain  the 
data  for  this  study  are  reviewed  next.  Magnetic 
field  waveforms  were  obtained  from  a single  wide- 
band sensor  mounted  on  the  centerline  of  the  up- 


per fuselage  of  the  WC-130  approximately  seven 
meters  from  the  nose  of  the  aircraft  (Bi  in  Fig- 
ure 1).  The  sensor  was  an  EG&G  model  CML-7  cyl- 
indrical Moebius  loop  B-Dot  sensor  (5,6)  that 
produces  an  output  voltage  proportional  to  the 
derivative  of  the  magnetic  field  over  a 38  MHz 
bandwidth.  A second  B-Dot  sensor  (Bo)  was  posi- 
tioned adjacent  to  the  first  and  oriented  to  re- 
spond to  magnetic  flux  parallel  to  the  wing  axis 
as  shown  in  Figure  1.  This  sensor  was  used  in 
conjunction  with  Bj  and  an  e.lectric  field  senior 
to  assist  in  locating  the  lightning  strikes.  The 
arrows  in  Figure  1 indicate  the  direction  of  in- 
creasing field  for  a positive  sensor  output  vol- 
tage. 

B-Dot  sensor  outputs  were  transmitted 
through  Merit  fiber  optic  data  links  to  instru- 
mentation inside  the  aircraft.  The  signals  were 
logarithmically  compressed,  digitized,  and  encod- 
ed on  magnetic  tape.  A partial  block  diagram  of 
tne  instrumentation  is  shown  in  Figure  2.  The 
digitizing  was  performed  with  a Micro  Pro  model 
PR7901  digital  transient  recorder  with  a 20  ns 
sample  period,  8 bit  amplitude  resolution  and  164 
ps  sample  window.  Forty  microseconds  of  pretrig- 
ger sampling  was  programmed  into  the  digitizer  so 
that  leader  pulses  preceeding  return  strokes 
could  be  recorded. 

An  electric  field  sensor  was  mounted  direct- 
ly ahead  of  the  B-Dot  sensors  on  the  upper  fuse- 
lage of  the  aircraft.  The  sensor  output  was 
integrated  electronically  and  recorded  on  two 
wideband  FM  channels  of  a Honeywell  model  101 
analog  tape  recorder.  Recorded  bandwidth  of  the 
electric  field  was  0.1  Hz  to  500  kHz.  The  analog 
recorder  also  recorded  a logic  pulse  indicating 
the  time  of  digital  transient  recorder  triggering. 
Trigger  pulse  location  within  the  overall  elec- 
tric field  waveform  was  the  first  criteria  used 
to  categorize  the  digitized  data  (Figure  3).  Re- 
finement of  the  categorization  was  obtained 
through  examination  of  the  fine  structure  of  the 
digitized  waveform. 

AIRCRAFT  EFFECTS  - We  have  simplified  the 
interpretation  of  the  magnetic  field  measurements 
to  a certain  degree  by  careful  positioning  of  the 
B-Dot  sensor.  By  positioning  the  sensor  in  the 
plane  defined  by  the  wings  and  upper  fuselage 
surface,  and  orienting  it  to  respond  to  magnetic 
flux  parallel  to  the  fuselage  axis  and  perpendic- 
ular to  the  wing  structure,  we  have  made  it  in- 
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sensitive  to  magnetic  fields  and  resonances  asso- 
ciated with  axial  wing  and  fuselage  current  flow. 
Some  coupling  of  resonant  responses  is  evidenced 
in  the  time  domain  waveforms  due  to  non-ideal 
symmetry  (8). 

Field  enhancement  due  to  the  aircraft  struc- 
ture can  be  determined  in  the  manner  described  by 
Taylor  (9,10).  The  aircraft  fuselage  is  modeled 
as  a cylinder  with  the  axial  component  of  the 
magnetic  field  at  the  surface  given  by 

Hj  » H^nc  (l-j2kacosfl) 

where  is  the  incident  magnetic  field,  a is 

the  fuselage  radius,  k«2ir/X  and  6 is  measured  up- 
ward from  the  direction  of  incidence  along  the 
fuselage  circumference.  This  equation  is  valid 
fcr  broadside  incidence  where  (ka)2<<l.  For  most 
of  our  data,  the  aircraft  is  positioned  so  that 
the  cosine  function  in  the  second  term  goes  to 
zero.  In  the  worst  case,  0 would  be  45°  and  the 
second  term  would  result  in  an  enhancement  of 
high  frequencies.  This  enhancement  would  be  be- 
low 3 dB  at  20  MHz,  however,  and  less  at  lower 
frequencies.  Therefore,  the  resultant  magnetic 
field  parallel  to  the  fuselage  axis  at  the  Bi 
sensor  location  is  approximately  equal  to  the  in- 
cident magnetic  field  at  this  point. 

RANGE  DETERMINATION  - Distances  to  thunder- 
storms were  determined  from  radar  patterns  ob- 
tained during  flight.  The  radar  data  were  re- 
corded on  digital  tape  and  processed  later  by 
computer  for  analysis.  Lightning  data  were  se- 
lected from  storm  cells  that  could  be  clearly 
identified.  When  several  storm- cells  were  active 
around  the  aircraft,  ambiguities  in  locations 
were  resolved  by  examining  the  polarities  and  re- 
lative amplitudes  of  the  B-Dct  signals.  Since 
the  precise  location  of  the  return  stroke  within 
storm  cells  was  unknown,  lightning  location  ac- 
curacy deteriorated  with  decreasing  thunderstorm 
range. 

DATA  PROCESSING 

The  digitized  lightning  data  was  scaled, 
anti-logged  and  numerically  integrated  to  display 
the  field  waveshape.  A typical  overall  data  win- 
dow and  the  integration  are  shown  in  Figure  3.  A 
10.24  us  (512  sample  points)  time  interval  con- 
taining the  data  of  Interest  was  selected  from 
the  overall  waveform  for  processing.  For  first 
return  strokes,  the  interval  was  chosen  to  in- 
clude the  alow  front  of  the  field  waveform  but  to 
exclude  stepped  leader  pulses.  Consequently,  the 
derivative  impulse  corresponding  to  the  fast 
front  of  the  return  stroke  was  placed  about  2 to 

4 us  after  the  beginning  of  the  interval.  Figure 

5 shows  a first  return  stroke  derivative  waveform 
and  the  numerical  integration  recorded  at  9 kilo- 
meters. Subsequent  return  stroke  waveforms  were 
positioned  approximately  1 to  2 us  after  the  be- 
ginning of  the  10.24  us  interval.  For  stepped 
leader  pulses,  4 u>  of  data  containing  the  pulse 
was  positioned  at  the  beginning  of  the  interval 
and  the  remaining  6.24  us  set  to  zero.  Subse- 


quent stroke  and  stepped  leader  waveforms  are 
shown  in  Figures  6 and  7,  respectively. 

A fast  Fourier  transform  of  the  data  was 
calculated  using  a rectangular  window  according 
to  the  following  equation: 

N^l  -J  (2itmn/N) 

X(mF)  - t£  x(nT)e 
n*0 

where  T is  the  sample  interval,  N is  the  number 
of  samples,  F = 1/NT,  m « 0,1,...,  N/2,  x(nT) 
is  the  sampled  magnetic  field  data,  and  X(mF)  is 
the  complex  spectra  value  (7).  With  512  samples 
and  a 20  ns  sample  period,  spectra  values  are 
calculated  at  discrete  frequency  increments  of 
approximately  97  kHz.  Since  field  derivatives 
were  recorded,  the  time  waveform  was  near  zero  at 
the  beginning  and  end  of  the  10.24  us  interval, 
minimizing  errors  introduced  by  abrupt  transi- 
tions at  the  boundaries  of  the  window.  Magni- 
tudes of  the  frequency  spectia  were  normalized  to 
a distance  of  50  km  assuming  a 1/r  relationship, 
where  r is  the  distance  from  the  thunderstorm  to 
the  aircraft,  and  corrected  for  angle  of  inci- 
dence to  the  sensor  by  dividing  by  the  sine  of 
the  angle  from  the  aircraft  nose  to  thunderstorm 
location.  Spectra  magnitudes  were  then  averaged 
at  each  frequency  and  standard  deviations  calcu- 
lated. The  derivative  response  of  the  sensor  was 
accounted  for  by  dividing  the  spectra  magnitudes 
by  radian  frequency  u>.  The  resulting  values  were 
expressed  in  decibels  using 

H (to)  (dB)  - 20  Lo£10  [|H(jm)|  ] 

where  the  horizontal  bar  denotes  mean  and 

the  vertical  bars  | | denote  magnitude.  Stand- 
ard deviations  (o(u>))  were  referenced  to  average 
values  and  expressed  in  decibels  using 

SD  (u)  (dB)  - 20  Log10  |H(jw)|  + o(u>)/|H(jm)  | j 
RESULTS 

The  data  were  recorded  during  several  storms 
occurring  on  25  August  and  26  August  1981  at  an 
altitude  of  5.2  km.  Ten  first  return  strokes  at 
a range  of  11  km  to  35  km  and  nine  subsequent  re- 
turn strokes  at  a range  of  7 km  to  28  km  were  se- 
lected for  processing.  The  last  stepped  leader 
pulse  preceeding  the  return  stroke  was  processed 
from  nine  of  the  first  stroke  waveforms.  Subse- 
quent strokes  showing  dart-stepped  leaders  were 
excluded  from  the  data  set. 

Average  spectra  values  for  first  return 
strokes,  subsequent  strokes  and  leader  pulses  ere 
plotted  in  Figures  8 through  10,  respectively.  A 
linear  interpolation  between  data  points  is  also 
plotted.  Spectra  magnitudes  st  selected  frequen- 
cies are  tabulated  in  Table  1 along  with  standard 
deviations. 

TRENDS  IN  SPECTRA  WITH  FREQUENCY  - Average 
spectre  values  st  100  kHz  are  -142  dB  for  both 
first  and  subsequent  return  strokes  and  20  dB 
lower  for  leader  pulses.  All  three  seta  of  data 
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decrease  approximately  inversely  with  frequency 
(f)  to  about  1 MHz  and  decrease  more  rapidly  at 
higher  frequencies.  Approximate  trends  in  aver- 
age spectra  amplitude  with  frequency  are  summar- 
ized in  Table  2. 

First  stroke  spectra  show  a decrease  propor- 
tional to  1 / f 2 from  800  kHz  to  5 MHz.  Beyond  5 
MHz  there  is  a slight  dip  in  the  spectra  and  then 
the  amplitudes  tend  toward  1/f  dependence. 

Subsequent  stroke  spectra  maintain  the  ini- 
tial 1/f  slope  from  100  kHz  out  to  1.2  MHz  and 
then  decrease  rapidly  until  3.5  MHz.  This  spec- 
tra also  regains  a 1/f  slope  by  20  MHz. 

Leader  pulses  show  less  relative  energy  at 
low  frequencies  with  a decrease  of  16  dE  from 
100  kHz  to  1 MHz  rather  than  20  dB  as  in  the  case 
of  the  return  strokes.  Beyond  1 MHz  the  spectra 
decreases  proportional  to  1/f. 

Evidence  of  aircraft  resonances  is  visible 
in  the  graphs  in  the  3 to  5 MHz  region. 

STANDARD  DEVIATIONS  - Standard  deviations 
are  greatest  for  leader  pulse  spectra  and  smal- 
lest for  subsequent  stroke  spectra  as  indicated 
in  Table  1.  The  larger  leader  pulse  standard 
deviations  are  indicative  of  the  more  random 
structure  of  the  leader  waveshapes  as  well  as  a 
lower  signal  to  noise  ratio.  Since  leader  pulse 
data  were  extracted  from  the  return  stroke  wave- 
forms, the  signal  to  noise  ratio  would  be  a fac- 
tor of  two  to  five  lower  for  leader  pulses  than 
for  first  return  strokes  (for  derivative  magnetic 
field  data).  Subsequent  stroke  standard  devia- 
tions are  iow,  particularly  at  lower  frequencies, 
suggesting  uniformity  in  different  return  stroke 
waveshapes. 

COMPARISON  TO  PUBLISHED  SPECTRA  - Frequency 
spectra  derived  from  broadband  ground-based  elec- 
tric field  measurements  have  been  obtained  by 
Serhan  (11)  and  more  recently  by  Weidman  (12). 
Serhan  calculated  spectra  from  2 kHz  to  700  kHz 
for  measurements  recorded  from  1.5  km  to  200  km. 
Weidman  obtained  spectra  of  electric  fields  from 
100  kHz  to  20  MHz  for  measurements  at  30  to  50  km 
where  propagation  was  over  saltwater  so  that  the 
effects  of  ground  wave  attenuation  were  minimized. 

Weidman' s data  for  first  return  strokes  show 
a spectral  amplitude  of  about  -95  dB  at  100  kHz, 
a slope  of  1/f  to  2 MHz,  1/f2  to  10  MHz,  and  1/f® 
beyond  10  MHz.  These  trends  are  similar  to  the 
WC-130  data  with  two  exceptions  — the  transition 
from  slopes  cf  1/f  to  1/f2  occurs  at  800  kHz  ra- 
ther than  2 MHz;  and  the  data  shows  no  rapid  de- 
crease in  spectral  amplitude  beyond  10  MHz  for 
the  aircraft  data. 

A comparison  in  spectral  amplitudes  between 
the  air  and  ground-based  data  is  difficult  for 
several  reasons.  For  Identical  waveshapes,  spec- 
* tra  amplitudes  should  be  proportional  to  waveform 

amplitudes.  With  different  frequency  characteris- 
tics, however,  the  proportionality  beconss  a 
function  of  frequency.  Also,  aircreft  waveforms 
are  obtained  at  varying  range  and  elevation  from 
the  return  stroke  which  precludes  a direct  com- 
parison with  ground-based  data  at  a fixed  posi- 
tion. Examination  of  the  first  stroke  data 
shows  a ratio  of  112  to  1120  for  electric  field 


spectra  on  the  ground  versus  magnetic  field  on 
the  aircraft  over  a 500  kHz  to  5 MHz  frequency 
range.  This  bounds  the  impedance  of  free  space 
(3770)  that  would  be  expected  for  simultaneous, 
distant,  ground-based  electric  and  magnetic  field 
data  (13). 

SUMMARY 

We  have  calculated  the  frequency  spectra  of 
broadband  lightning  magnetic  field  measurements 
recorded  on  a WC-130  aircraft  in  flight.  The 
data  were  recorded  at  a range  of  7 to  35  km. 
Spectra  values  were  calculated  from  100  kHz  to 
20  MHz  for  stepped  leader  pulses,  first  return 
strokes  and  subsequent  return  strokes.  Spectra 
amplitude  trends  with  frequency  were  extracted 
and  compared  with  trends  of  published  ground 
based  spectra.  From  the  data  we  can  infer  the 
following  lightning  characteristics: 

1.  First  and  subsequent  return  strokes 
have  similar  energy  levels  in  the  100  kHz  to  20 
MHz  frequency  range. 

2.  The  initial  transition  point  from 
slopes  of  1/f  to  steeper  slopes,  occurring  st  800 
kHs  for  first  strokes  and  1.2  MHz  for  subsequent 
strokes,  is  probably  indicative  of  faster  subse- 
quent stroke  risetimes. 

3.  Leader  spectra  lie  14  to  20  dB  be- 
low first  stroke  spectra  at  low  frequencies  sug- 
gesting a factor  of  5 to  10  difference  in  magnet- 
ic field  amplitude. 

4.  Subsequent  stroke  spectra  show  a 
relatively  small  standard  deviation  suggesting  a 
uniformity  of  waveshape  among  independent  return 
strokes.  Although  some  of  these  characteristics 
could  be  observed  from  time  domain  waveforms,  one 
significant  feature  that  is  extracted  by  Fourier 
analysis  is  the  presence  of  high  frequency  energy 
beyond  10  MHz  that  is  not  observed  with  ground- 
based  measurements.  From  our  analysis,  this  does 
not  appear  to  be  characteristic  of  the  recording 
system  or  airframe  enhancement  but  rather  an  in- 
herent feature  of  the  lightning  fields. 
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Table  1 - Average  Spectra  Values  and  Standard  Deviations  in  Decibels 
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Frequency 

First 

Mean 

Strokes 

SD 

Subsequent 

Mean 

Strokes 

SD 

Leaders 
Mean  SD 

100  kHz 

-192 

2.4 

-142 

2.1 

-162 

4.2 

200  kHz 

-152 

4.5 

-146 

2.2 

-166 

3.5 

500  kHz 

-158 

3.7 

-156 

2.8 

-172 

3.1 

1MHz 

-170 

3.1 

-164 

2.5 

-179 

4.0 

2MHz 

-177 

3.4 

-177 

3.0 

-191 

4.3 

5 K z 

-196 

4.2 

-196 

4.0 

-205 

4.2 

10  MHz 

-210 

3.5 

-209 

3.3 

-216 

3.9 

15  MHz 

-213 

2.9 

-215 

5.2 

-220 

3.6 

20  MHz 

-218 

4.6 

-221 

4.6 

-225 

4.7 

Table  2 - Variations  of  Spectra  Amplitude  with  Frequency 


First  Return  Strokes 

Subsequent  Return  Strokes 

Leader  Pulses 

Frequency  Range 

Trend 

Frequency  Range 

Trend 

Frequency  Range  Trend 

0.1  - 0.8  MHz 
0.8  - 5.0  MHz 
5.0  - 20  MHz 

1/f 

l/f2 

1/f 

0.1  - 1.2  MHz 
1.2  - 3.5  MHz 
10  - 20  MHz 

1/f, 

1/f3 

1/f 

0.1-1.0MHz  -16dB/decade 
1.0  - 20  MHz  1/f 
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Fig.  4 - Overall  graph  of  digital  magnetic  field 
erivative  data  (H,  top  trace)  and  the  numerical 
integration  (bottom  trace) 
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Fig.  5 - Ten  microsecond  data  window  for  first 
return  stroke  waveform  at  9.5  km  with  numerical 
integration 
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Fig.  9 - Graph  showing  average  spectra  amplitude 
in  dicibels  versus  log  frequency  for  subsequent 
return  strokes 
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Fig.  7 - Ten  microsecond  data  window  for  leader 
pulse  waveform  at  approximately  11  km  with  numer- 
ical integration 
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Fig.  8 - Graph  showing  average  spectra  amplitude  Fig.  10  - Graph  showing  average  spectra  amplitude 
in  decibels  versus  log  frequency  for  first  return  in  decibels  versus  log  frequency  for  leader 
strokes  pulses 
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